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Motivated by results implying that the constituents of dark matter (DM) might be coUisional, we consider a cosmological 
(toy-) model, in which the DM itself possesses some sort of thermodynamic properties. In this case, not only can the 
matter content of the Universe (the baryonic component, which is tightly gravitationally-bounded to the dark one, also 
being included) be treated as a classical gravitating fluid of positive pressure, but, together with all its other physical 
characteristics, the energy of this fluid's internal motions should be taken into account as a source of the universal 
gravitational field. In principle, this form of energy can compensate for the extra (dark) energy, needed to compromise 
spatial flatness, while the post-recombination Universe remains ever-decelerating. What is more interesting, is that, 
at the same time (i.e., in the context of the coUisional- DM approach), the theoretical curve representing the distance 
modulus as a function of the cosmological redshift, ^i{z), fits the Hubble diagram of a multi-used sample of supernova 
la events quite accurately. A cosmological model filled with coUisional DM could accommodate the majority of the 
currently-available observational data (including, also, those from baryon acoustic oscillations), without the need for 
either any dark energy (DE) or the cosmological constant. However, as we demonstrate, this is not the case for someone 
who, although living in a Universe filled with self-interacting DM, insists on adopting the traditional, coUisionless-DM 
approach. From the point of view of this observer, the cosmologically-distant light-emitting sources seem to lie farther 
(i.e., they appear to be dimmer) than expected, while the Universe appears to be either accelerating or decelerating, 
depending on the value of the cosmological redshift. This picture, which, nowadays, represents the common perception in 
observational cosmology, acquires a more conventional interpretation within the context of the coUisional-DM approach. 
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1. Introduction 

The beginning of the 21** century was one of the most ex- 
citing epochs for cosmology as a science. According to ob- 
servational data on the temperature variations in the cos- 
mic microwave background (CMB) that became in public 
at that epoch (de Bernardis et al. 2000; Jaffe et al. 2001; 
Padin et al. 2001; Stompor et al. 2001; Netterfield et al. 
2002), now, we are quite confident that the Universe can be 
adequately described by a spatially flat Robertson- Walker 
(RW) cosmological model. 

As a consequence, the total energy density, e, of the 
Universe matter-energy content, in units of the energy den- 
sity Ec = PcC^ (equivalent to the critical rest-mass density, 

Pc = 1^ , where Hq is the Hubble parameter at the present 
epoch, c is the velocity of light, and G is Newton's univer- 
sal constant of gravitation), should be very close to unity, 
O = — ~ 1, i.e., much larger than the measured quantity, 



n 



M 



0.3 (Komatsu at al. 2009). 



At the same time, high-precision distance measure- 
ments, performed with the aid of the supernovae la (SNe 
la) standard candles, indicated that, in any cosmological 
model with vanishing cosmological constant. A, the far-off 
light-emitting sources appear to be dimmer than expected 
(Riess et al. 1998; Perlmutter et al. 1999). 
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The observational data then seemed to favour a 
Universe of collisionless content (i.e., filled with matter in 
the form of dust) and A ^ 0, in which £7^/ — 0.3 and 

17a = 1^ ~ 0.7 (Riess et al. 2001, 2004). Since a non- 

vanishing cosmological constant (necessarily) involves a re- 
pulsive (gravitational) force (see, e.g., Sahni 2004), the ap- 
parent dimming of the distant light-emitting sources was 
attributed to a relatively recent phase of accelerated ex- 
pansion (see, e.g., Linder 2008). 

The onset of the dimming of the cosmologically distant 
indicators used (hence, the associated transition from ac- 
celeration to deceleration), takes place at a relatively low 
value of the cosmological redshift, z, the so-called transition 
redshift, zt, which, nowadays, is being (observationally) set 
at Zt = 0.46 ± 0.13 (Riess et al. 2004, 2007). In this case, 
the cosmological constant can be determined observation- 

ally, since, on theoretical grounds, zt = (2j^ j — 1 (see, 

e.g., Perivolaropoulos 2007, Eq. (21)), i.e., the transition 
redshift depends on the value of A. The particle-physics 
vacuum does contribute an effective cosmological constant, 
which could serve (also) as compensation to the extra en- 
ergy needed to flatten the Universe (Sahni & Starobinsky 
2000). Unfortunately, the energy-density attributed to such 
a source is 10^^'^ times larger than what is observed (see, 
e.g., Padmanabhan 2003; Sahni 2004). 
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Hence, it became evident that, for the above-mentioned 
observational results to be reconciled vifithin a unified theo- 
retical framework, a different approach (i.e., other than the 
cosmological constant) was needed. 

In this context, theorists have focused on two main ar- 
eas: (i) The introduction of an exotic, negative-pressure 
fluid, the dark energy (Caldwell et al. 1998), occasionally re- 
ferred to as quintessence (Carroll 1998), and (ii) the con- 
sideration of alternative-gravity theories, such as the scalar- 
tensor theories (Esposito-Farese & Polarski 2001) and the 
/(i?)-gravity (Capozzicllo et al. 2003), together with sev- 
eral braneworld scenarios, including DGP-gravity (Dvali et 
al. 2000) and the landscape scenario (Bousso & Polchinski 
2000). 

Other physically-motivated models, predicting an ac- 
celerated expansion, have also appeared in the literature, 
involving holographic gravity (Cohen et al. 1999; Li 2004; 
Pavon & Zimdahl 2005), Chaplygin gas (Kamenshchik ct al. 
2001; Bean & Dore 2003; Sen & Scherrer 2005), Cardassian 
cosmology (Freese & Lewis 2002; Wang et al. 2003), theo- 
ries of compactified internal dimensions (Perivolaropoulos 
2003), and mass- varying neutrinos (Fardon et al. 2004; 
Peccei 2005). 

However, most of these attempts were inhibited by the 
so-called coincidence problem, i.e., the need to explain why 
the Universe transferred from deceleration to acceleration 
so recently (see, e.g., Perivolaropoulos 2007). 

It has been more than a decade since the first ob- 
servations, which provoked the aforementioned scientific 
(r)cvolution, became in public, and, still, no undisputed 
theoretical framework has been developed to accommo- 
date them (for a detailed review see, e.g., Caldwell & 
Kamionkowski 2009). In the meantime, the need for an ex- 
tra (dark) energy component has been confirmed by other 
observational methods, including galaxy clusters dynamics 
(Allen et al. 2004), the integrated Sachs- Wolfe (ISW) effect 
(Boughn & Crittenden 2004), and baryon acoustic oscilla- 
tions (BAG) (Eisenstein et al. 2005; Percival et al. 2010). 
In spite of the wealth of references, previous studies have 
been far from exhaustive (see, e.g., Albrecht et al. 2006; 
Peacock et al. 2006) and, perhaps, we should keep our op- 
tions open, also, to more conventional interpretations (see, 
e.g., Buchert 2000, 2001; Kolb et al. 2006; Celerier 2007; 
Ellis 2009). 

On the other hand, much evidence has already been 
accumulated in support of a (non-baryonic) DM compo- 
nent in the Universe matter-content (see, e.g., Tegmark 
et al. 2006; Spergel et al. 2007). Among other pieces of 
evidence, this support includes flattened galactic rotation 
curves (Begeman et al. 1991; Borriello & Salucci 2001), 
the weak gravitational Icnsing of distant galaxies by (some 
dark) foreground structure (Hoekstra ct al. 2002), and the 
weak modulation of strong lensing around individual mas- 
sive elliptical galaxies (Moustakas & Metcalf 2003). On the 
scale of galaxy clusters, observations (of radial velocities, 
weak lensing, and X-ray emission) indicate a total mass 
density almost ten times higher than the corresponding 
density in baryons (Bahcall & Fan 1998; Kashlinsky 1998; 
Tyson et al. 1998). On cosmological scales, the anisotropics 
observed in the CMB have led to an estimate of the total- 
mass density of the order of J^m^^ = 0.1358±0.0037, where 
h is the Hubble parameter in units of 100 km/ sec per Mpc 
(Komatsu et al. 2009). In contrast, measurements of the 
light-chemicals' abundances (Olive et al. 2000) have led to 



an estimate of the baryonic-mass density of the order of 
nsh'^ = 0.02273 ± 0.00062. The combination of these re- 
sults suggests that, more than 85% (by mass) of the matter 
in the Universe consists of non-luminous and non-baryonic 
material. 

Although we do not know for certain how the DM 
came to be formed, a sizeable relic abundance of weakly- 
interacting massive particles (WIMPs) is generally ex- 
pected to have been produced as a by-product of the 
Universe's hot youth (see, e.g., Kolb & Turner 1990, p. 369). 
These particles decouple from radiation much earlier than 
pure-baryonic matter does. Hence, very soon after recombi- 
nation (tfi), the baryons fall into deep potential wells of the 
already evolved DM-perturbations and become bounded to 
them, i.e., for t > Ir, there arc no freely-floating baryons 
around (Olive 2003). 

Among the various candidates for DM constituents, the 
thermal WIMPs remain one of the most attractive. They 
appear, generically, in theories of weak-scale physics be- 
yond the standard model, while giving the appropriate relic 
abundance (Srednicki et al. 1988; Gondolo & Gelmini 1991). 
These particles are also helpful to consider in terms of di- 
rect and indirect detection of DM (see, e.g., Jungman et al. 
1996; Bertone et al. 2005), because they must have some 
connection to standard- model particles (see, e.g., Hooper 
2009). In addition to debating their precise nature, the sci- 
entific community used to argue that the WIMPs should 
be coUisionless. 

However, many results from high-energy particle de- 
tectors, such as the ATIC (Chang et al. 2008) and 
PAMELA (Adriani et al. 2009), combined with data from 
the Wilkinson microwave anisotropy probe (WMAP) sur- 
vey (Hooper et al. 2007), have revealed an unusually high 
electron - positron production in the Universe, much more 
than anticipated by SNc explosions or cosmic-ray collisions. 
These results have led many scientists to argue that among 
the best candidate sources of these high-energy events are 
the annihilations of WIMPs (see, e.g., Barger et al. 2008; 
Bergstrom et al. 2008; CireUi & Strumia 2008; Regis & Ullio 
2008; Baushev 2009; Chohs et al. 2009a, 2009b; Fornasa et 
al. 2009; Fox & Poppitz 2009; Kane et al. 2009; Zurek 2009, 
for an extensive, though incomplete list), i.e., that the DM 
constituents can be slightly collisional (see, e.g., Spergel & 
Steinhardt 2000; Arkani-Hamed et al. 2009; CireUi et al. 
2009; Cohen & Zurek 2010), although, some studies dis- 
agree with this interpretation (see, e.g., Feng et al. 2010). 

A cosmological model filled with self-interacting DM 
could be a relatively inexpensive solution to the DE prob- 
lem, and, several ways of accommodating both the DM 
and the DE into a unified theoretical framework have 
been considered (see, e.g., Zimdahl et al. 2001; Bilic et al. 
2002; Balakin et al. 2003; Scherrer 2004; Lima et al. 2008; 
Basilakos & Plionis 2009, 2010; Dutta & Scherrer 2010). 
In this context, wc suggest that, phenomenologically, the 
self-interacting DM could attribute to the Universe matter- 
content some sort of fluid-like properties, and (so) lead to 
a conventional approach to the DE concept. 

The main outstanding problem of the current cosmolog- 
ical picture is that the Universe must contain an amount of 
energy that is considerably higher than the equivalent of the 
total rest-mass of its matter content. However, if the DM 
constituents collided with each other frequently enough, en- 
abling their (kinetic) energy to be re-distributed, i.e., if the 
DM itself possessed some sort of thermodynamic properties, 
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a conventional extra-energy component might be present in 
the Universe, given by the energy of the internal motions 
of the coUisional-DM fluid. 

On this basis, it is worth examining the evolution and 
the dynamical characteristics of a cosmological model (not 
necessarily reflecting our own Universe), in which (in princi- 
ple) there is no DE at all. The matter-energy content of this 
model consists solely of two components, i.e., the DM (dom- 
inant) and the baryonic one (subdominant), both having 
the abundances attributed to them by analyses of the five- 
year survey of the WMAP (Dunkley et al. 2008; Komatsu 
et al. 2009). Accordingly, we demonstrate that these two 
components are (by themselves) sufficient (i) to repro- 
duce the result that (today) ft = 1, (ii) to account for 
the observed dimming of the distant light-emitting sources, 
and (iii) to explain the apparent accelerated expansion 
of the Universe. All the above provided that, macroscopi- 
cally, these two constituents (basically the dark one) form 
a gravitating fiuid with a tliermodynamical content. In this 
case, together with all the other physical characteristics, 
the energy of this fluid's internal motions should (also) be 
taken into account as a source of the universal gravitational 
field. Although speculative, the idea that the extra (dark) 
energy needed to flatten the Universe could be attributed 
to the internal motions of a collisional-DM fluid is (at least) 
intriguing. 

This paper is organized as follows. In Sect. 2, we 
explore the dynamical characteristics of a cosmological 
model driven by an ideal fluid, consisting (mainly) of 
(thermodynamically-involved) DM, with positive pressure, 
the volume elements of which perform adiabatic flows. 
Accordingly, after deriving the corresponding scale factor, 
we determine the functional form of several parameters of 
cosmological significance, each one depending on the cos- 
mological rcdshift, such as the luminosity distance and 
the distance modulus of the cosmologically-distant light- 
emitting sources, together with the Hubble and the decel- 
eration parameters, which characterize the cosmic expan- 
sion. The corresponding results suggest that, in the con- 
text of the collisional-DM treatment, the extra (dark) en- 
ergy (needed to compromise spatial flatness) can be com- 
pensated by the energy of the internal motions of this 
fluid, while, the post-recombination Universe remains ever- 
decelerating. However, as we demonstrate in Sect. 3, this 
is not the case for someone who (although living in a 
Universe filled with coUisional DM) insists on adopting the 
traditional (collisionless-DM) approach. From the point of 
view of this observer, besides the need for an extra-energy 
component (for confronting the CMB-based observational 
results), the cosmologically-distant light-emitting sources 
seem to lie farther (i.e., they appear to be dimmer) than ex- 
pected, while the Universe appears to be either accelerating 
or decelerating, depending on the value of the cosmological 
rcdshift. Finally, we conclude in Sect. 4. 



2. A Universe filled with collisional dark matter 

It is generally accepted that the study of the CMB 
has proven to be a powerful tool in exploring the post- 
recombination Universe. According to the various CMB- 
oriented observational data, the Universe has emerged out 



of the radiation epoch as a spatially-flat RW model (see, 
e.g., de Bernardis et al. 2000) 



ds^ = S^{t]) [c^dif - {dx^ + dy^ + dz^)] 



(1) 



where 77 is the conformal time and S'(r/) is the scale factor. 
As a consequence, the value of the Hubble parameter at the 
present epoch is, by deflnition, given by 



SttG 



(2) 



(see, e.g.. Peacock 1999, p. 77). The evolution of this model 
depends on the nature of the source that drives the univer- 
sal gravitational field, i.e., its matter-energy content. 

Along the lines of the collisional-DM approach, in spec- 
ifying the Universe matter-energy content, wc assume that, 
in principle, there is no DE at all. Instead, we admit that 
the DM, together with the small, baryonic "contamination" 
(the latter is too tightly gravitationally-bounded to the for- 
mer), possess fluid-like properties. In this sense, the colli- 
sions of the WIMPs maintain a tight coupling between them 
and their energy can be re-distributed, i.e., the DM itself 
also possesses some sort of thermodynamical content. In 
this case, the evolution of the post-recombination Universe 
is no longer driven by dust, but, by a gravitating fluid of 
positive pressure, p, satisfying the equation of state 



p — wpc 



(3) 



where p is the rest- mass density (the part, equivalent to the 
energy density pc^ , that remains unaffected by the internal 
motions of the cosmic fluid) and 0<w=(^) <lisa 
dimensionless constant, which measures the square of the 
speed of sound, c^, in units of c^. Now, the fundamental 
units of the Universe matter-content are the volume ele- 
ments of the collisional-DM fluid (elements of fluid, each 
one consisting always of the same particles). 

The motions of the volume elements in the interior of a 
continuous medium are governed by the equations 



T^r. = 



(4) 



where Greek indices refer to the four-dimensional space- 
time (in connection, Latin indices refer to the three- 
dimensional spatial slices), the semicolon denotes covariant 
derivative, and T'^^ is the energy-momentum tensor of the 
Universe matter-content, i.e., basically (but not solely), of 
the collisional-DM fluid. 

Conflning ourselves to the particular case of a perfect 
fluid, T^"^ takes on the standard form 



Tf"" = {e+ p)uf'u'' - pg^"" , 



(5) 



where = dx'^/ds is the four- velocity {u^u^ = 1) at the 
position of a fluid's volume element, g'^^ are the contravari- 
ant components of the Universe metric tensor, and e is this 
fluid's total-energy density. In an (ideal) equilibrium state, 
i.e., in the absence of shear, viscocity and heat conductivity, 
e is decomposed to 



e = pc + pH 



(6) 



(for a detailed analysis see, e.g., Fock 1959, pp. 81 - 83 and 
91 - 94), where H is the potential energy per unit rest-mass, 
associated with the infinitesimal deformations (expansions 
or/and compressions) of the fluid. Upon consideration of 
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adiabatic processes, 11 coincides with the energy of this 
fluid's internal motions (per unit rest-mass), thus defining 
pll as the corresponding energy density, associated with the 
thermodynamical content of the (collisional) DM. 

Along these lines, the equations represented by Eq. (4) 
are the hydrodynamic flows of the volume elements in the 
interior of a perfect-fluid source 



1 



P 



(7) 



where the comma denotes a partial derivative and (5^ is the 
Kronecker symbol. The equations given in Eq. (7) can be 
cast in the more convenient form 



ds 



1 



(8) 



where V^^^ are the Christoffel symbols corresponding to the 



Universe metric tensor, g^^^, and h'^^ 
projection operator. 

However, in a maximally symmetric cosmological setup, 
there is no real difference between hydrodynamic flows and 
the ballistic motions along the (non-intersecting) geodesic 
trajectories a;^ = constant, since, the equations of motion 
given by Eq. (8) are trivially satisfied by the fluid's volume 



elements with ' 



(1, 0, 0, 0). In other words, in comoving 



coordinates, both sides of Eq. (8) vanish identically, result- 
ing in the geodesic equations 







and the equations of the ballistic trajectories 

/i"^P,A = . 



(9) 



(10) 



Notice however that, in this case, Eq. (10) is satisfied only 
because of the form of the metric tensor in Eq. (1), i.e., even 
if \ does not vanish. Hence, in comoving coordinates, the 
geodesic motions and the hydrodynamic flows of the cos- 
mological model given by Eq. (1) are, practically, indistin- 
guishable. Therefore, a comoving observer of the cosmic 
expansion also traces the hydrodynamic flow of the homo- 
geneous cosmic fluid and the Weyl's postulate is valid (see, 
e.g., Narlikar 1983, p. 91). 

As a consequence, the dynamical evolution of the model 
given by Eq. (1) is governed by the Friedmann equation 
(with A = 0) of the classical Friedmann-Robertson- Walker 
(FRW) cosmology 



SttG 



where 



H = 



(11) 



(12) 



is the Hubble parameter as a function of the scale factor, 
and the prime denotes differentiation with respect to 77. 

Nevertheless, inherently, there is an essential difference 
between our model and the rest of the classical FRW cos- 
mologies. In our case, the basic matter constituents (al- 
though they may resemble test particles receding from each 
other) are the volume elements of a collisional-DM fluid, 
i.e., they possess some sort of internal structure, hence ther- 
modynamical content. Therefore, the functional form of s 



in Eq. (11) is no longer given by pc^ alone, but by Eq. (6) 
(see also Narlikar 1983, pp. 61, 62). 

In this model, the first law of thermodynamics for adi- 
abatic flows, given by. 



dn+pd = 

(see, e.g., Chandrasekhar 1965), results in 

n = Ho + wc^ In ( — 
\Po 



(13) 



(14) 



where the constants po and Hq are assumed to denote the 
corresponding present-time values. Accordingly, the total- 
energy density of the Universe matter-energy content is 
written in the form 



u'^u^ is the ^ ^'^ 



1 



no^ , fp 

— +wln[ — 



(15) 



On the other hand, for every value of if, the conservation 
law T^^ ~ 0, in terms of the metric tensor of Eq. (1), yields 



£' + 3^(e + p) = 0, 



(16) 



which, upon consideration of Eqs. (3) and (6), results in 



, ^0 

P = Po\ Y 



(17) 



where Sq is the value of S'(r/) at the present epoch. Equation 
(17) represents the conservation of the total mass in a cos- 
mological model in which matter dominates, i.e., for every 
T] within the post-recombination epoch (see, e.g., Tsagas et 
al. 2008). 

With the aid of Eqs. (15) and (17), Eq. (11) is written 
in the form 



SttG 



-Po 



1 



Ho 

^2 



3w In 



Now, combining Eqs. (2) and (18), we obtain 



H 



Ho 

r.2 



?)W In 



5, 



(18) 



(19) 



At the present epoch, where S = Sq and H = Hq, we have 



n 



Ho 



1, 



(20) 



from which, the present-time value of the internal energy 
per unit rest- mass, Hq, emerges as 



Hn 



1 



1 



(21) 



Since flM < 1, Eq. (21) suggests that, at the present epoch, 
the energy density, PqHq, of the internal motions of a grav- 
itating perfect fluid (consisting, mainly, of collisional DM 
and a small baryonic contamination) dominates over the 
corresponding rest-mass quantity, i.e.. 



Polio = ( ^ - 1 ) poc^ > Poc^ 



(22) 
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But, what is more important, is that, at the same time, the 
combination of Eqs. (6) and (21) results in the following 
value of the total-energy density parameter 



PcC^ 



PcC^ 



= O 



M 



M- 



= 1 



(23) 



In view of Eq. (23), the extra (dark) energy, needed to flat- 
ten the Universe, can be provided by the energy of the in- 
ternal motions of a thcrmodynamically-involved-DM fluid. 

On the other hand, upon consideration of Eq. (21), Eq. 
(19) is written in the form 



H 



1 + SwflM 111 ( ^ 



(24) 



Equation (24) can be solved, explicitly, in terms of the error 
function (see Appendix A). However, it can become partic- 
ularly transparent (and useful) if we take into account that, 
since < w < 1 and CIm — 0.3, the combination wfl]\i can 
be quite small, i.e., w^Im ^ 1. In this case, we may take 
the natural logarithm on both sides of Eq. (24) , to obtain 



In 



H 



= In 



In 



1 + wCIm 111 



(25) 



Within the post-recombination era, ^ < 1090, hence 
In(^)'^ < 21. Therefore, as long as wJIm ^ 1, we have 



In 



1 4- wVIm In ( ^ 



wQAf In ( ^ 



so that, for terms linear in wflM, Eq- (24) results in 



(26) 



(27) 



In this case, using Eq. (12), we can solve Eq. (27), to de- 
termine the scale factor of the coUisional-DM model (1), as 
follows 



So 



(28) 



where we have defined the present-time value, rjQ, of the 
conformal time, t], as 



(1 + 3wnM)HoSo 



(29) 



For w ^ 0, Eq. (28) is the natural generalization of 
the corresponding collisionless-DM model, the well-known 
Einstein-dc Sitter (EdS) Universe (5* ^ 77^) (see, e.g.. 
Peacock 1999, pp. 77, 83 and 142 - 144). 

Eventually, in the collisional-DM model given by Eq. 
(1), the cosmological redshift parameter is defined as 



thus Eq. (27) is written in the form 
= i?o(l + 2:)^^^+"'""^ . 



(30) 



(31) 



We note the striking functional similarity between Eq. (31) 
and the corresponding result for a dark-energy fluid with 



equation of state in the form of Eq. (3) (cf. Eqs. (13) and 
(14) of Perivolaropoulos 2007). In our case, however, w > 
and; therefore, on the approach to z — 0, H{z) decreases 
monotonically. In other words, a cosmological model filled 
with collisional DM, necessarily, decelerates its expansion. 

This can be readily verified, by expressing the corre- 
sponding deceleration parameter, g, in terms of H and z. 



q{z) 



dH/dz 



(1 + ^) 



(32) 



(cf. Eq. (16) of Perivolaropoulos 2007), which, in view of 
Eq. (31), yields 



q{z) = -{i + SwriM) >o. 



(33) 



independently of z, even for w = 0. In other words, the 
model of a gravitating perfect-fluid source, as it stands, 
i.e., either pressureless (geodesic motions) or not (hydrody- 
namic flows), seems to be inappropriate for explaining the 
apparent accelerated expansion of the Universe. 

The actual reason is that it does not have to ac- 
count for any acceleration at all. As we demonstrate in 
the next Section, in a Universe filled with collisional (i.e., 
thermodynamically-involved) DM, the observed dimming 
of the distant light-emitting sources can be explained with- 
out the assumption of the accelerated expansion. 

3. Mistreating the dark matter as collisionless 

When the (unexpected) dimming of the SNc la standard 
candles was first discovered, the common perception about 
the cosmos, to the best of our knowledge, was that the 
DM constituents are collisionless, thus the various motions 
in the Universe were (necessarily) interpreted as geodesic 
motions of test particles receding from each other, i.e.. 



ds 



0. 



(34) 



Tilde variables are used, to distinguish the various quan- 
tities in Eq. (34) from the corresponding quantities used 
in Eq. (9), thus reflecting that the physical content of a 
collisionless-DM Universe (in which both the pressure and 
the energy of the internal motions are assumed to be negli- 
gible and, therefore, disregarded) is entirely different from 
that of the thcrmodynamically-involved-DM model (where 
p, H^O). 

In other words, the dynamical properties of a dust 
model are no longer described by 5^1/, i.e., Eq. (1), but 
rather in terms of another metric tensor, 5^1/, for which the 
corresponding (spatially-flat) line-element is written in the 
form 



dS^ = R^{t]) [c^dif - {dx^ + dy^ + dz^)] 



(35) 



Clearly, the evolution of this model is given in terms of the 
scale factor R{ri). From the point of view of an observer 
who (mis)treats the DM as collisionless, g^j, is the metric 
tensor upon which he/she should rely on, in interpreting 
observations. 

However, we recall that, for this observer, the accu- 
mulated evidence in favour of spatial flatness (necessarily) 
leads to the assumption of an extra (dark) energy com- 
ponent, in contrast to the collisional-DM case, where this 
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assumption would no longer be necessary. In the latter 
case, the appropriate candidate to provide the extra en- 
ergy needed to flatten the Universe is already included in 
the model (the energy of the internal motions). 

Furthermore, in the collisionless-DM scenario, every 
theoretical effort to interpret the (apparent) dimming of 
the SNe la standard candles, naturally, should also be based 
on g^u and the cosmologically relevant parameters arising 
from it. Accordingly, a possible explanation could be that, 
recently, the Universe accelerated its expansion (Riess et 
al. 1998; Perlmutter et al. 1999). This assumption, how- 
ever, attributes unnecessarily-exotic properties to the extra 
amount of energy needed to account for the spatial flat- 
ness (e.g., it should be repulsive in nature, i.e., of negative 
pressure, etc.). Therefore, we cannot help but wondering 
whether there is another (more conventional) explanation 
to be found (also) within the context of the coUisional-DM 
model. 

In what follows, we demonstrate that both the observed 
dimming of the distant light-emitting sources and the accel- 
erated expansion of the Universe could be only apparent, 
based on the misinterpretation of several cosmologically- 
relevant parameters, by someone who (although living in a 
Universe filled with coUisional DM) insists on adopting the 
traditional (collisionless-DM) approach. 

To explore this possibility, we note that the collisional- 
DM treatment of the Universe's matter content (in terms 
of which p ^ and the motions of its constituents are, 
in principle, hydrodynamic flows) can be related to the 
collisionless-DM approach where p = and the corre- 
sponding motions are, necessarily, geodesies, by means of 
a conformal transformation of the metric tensor (Kleidis & 
Spyrou 2000; Spyrou 2005; Spyrou & Tsagas 2004, 2010). 
Accordingly, from the original metric, g^,y, in terms of 
which the hydrodynamic flows have their well-known form 
given by Eq. (8), we can transfer the problem to a vir- 
tual metric, g^u, in terms of which the volume elements 
of the fluid move along geodesic-like trajectories, i.e., their 
velocity-vector obeys Eqs. (34) (in connection, see Synge 
1937; Lichnerowicz 1967, pp. 24 - 29 and 54-61; Carter 
1979). The appropriate (conformal) transformation for such 
a transition is 



(36) 



where, upon consideration of isentropic flows, the conformal 
factor F{x'^) takes on the functional form (Kleidis & Spyrou 
2000) 



Fix" 



C 



P 



C 



1 



1 n + 5 

c \ P 



(37) 



with C being an arbitrary (integration) constant. From Eq. 
(37), it becomes evident that F{x'^) is, essentially, the spe- 
ciflc enthalpy of the ideal fluid under consideration. 

Verozub (2008) extrapolated these results to include 
every Riemannian space-time and not just the metric at- 
tributed to a bounded, perfect-fluid source. In particular, 
he showed that the adiabatic hydrodynamic motion of an 
ideal- fluid element in a space-time with metric tensor g^^^, 
takes place along the geodesic lines of a Riemannian man- 
ifold with metric tensor given by the combination of Eqs. 
(36) and (37). 

With the aid of the technique developed by Kleidis & 
Spyrou (2000), we then determine the scale factor of the 



spatially-flat cosmological model in Eq. (35), i.e., the scale 
factor of the Universe as inferred by someone who, although 
living in a coUisional-DM Universe (where p, 11 > and 
^ 7^ 0), misinterprets the DM as coUisionless (p = 0). 

In principle, one can (always) use a (conformal) trans- 
formation to remove (from the rhs of Eqs. (8)) either 
the pressure gradient, which measures the response to 
non-gravitational forces, or the pressure itself (see, e.g., 
Lichnerowicz 1967, p. 26). Nevertheless, something like 
this (usually) comes with a price (see, e.g., Bruncton & 
Esposito-Farese 2007): The new metric, in terms of which 
the DM appears to be pressureless, is no longer a solution of 
general relativity (GR), but, rather, a solution to a modifled 
theory of gravity. The reason is that, in terms of this new 
metric (i.e., after the transformation (36) is applied), the 
action of the original gravitational fleld is also modified, 
acquiring extra terms in addition to the Einstein-Hilbert 
Lagrangian (see Appendix B). In other words, every effort 
to treat a collisional-DM fluid as pressureless, cannot be ac- 
complished in the context of GR. As a consequence, R{ri) 
is no longer a solution of the original Friedmann equation, 
given by Eq. (11). 

In view of Eq. (36), the scale factor of the Universe as 
it is inferred by a supporter of the coUisionlcss-DM sce- 
nario, R{r]), is related to the corresponding quantity of the 
collisional-DM model, S{'i]), as follows 



(38) 



where, by virtue of Eqs. (7), (17), and (30), F{x'^) is given, 
in terms of z, by 

F{z) = -^{1 + wnM [I + 3ln{l + z)]) . (39) 

In Eq. (39), the arbitrary integration constant, C, can be 
determined, by demanding that, in the (isobaric) pressure- 
less case, these two models should coincide, i.e., g^iy — 
In other words, for w = = p, R{r'i) = S{ri,w ~ 0), which 
represents the EdS model. Hence, F{w = 0) = 1. 

For p = constant = po, the first law of thermodynamics 
given in Eq. (13) yields 



n + — = constant 
P 



(40) 



which, in the particular case of dust, where po = Oj results 
in n = constant = Hq. Now, from Eq. (19) it becomes 
evident that, at the present epoch, 11 (= Hq) is (also) given 
by Eq. (21), even for w = Q. Accordingly, inserting Eq. (21) 
into Eq. (37), we find that, the condition F{w = 0) = 1 
leads to 



C = Q. 



M 



As a consequence, Eq. (39) results in 
F{z) = 1 + u;17a/ [1 + 3 ln(l + z)] . 



(41) 



(42) 



Using Eqs. (38) and (42), we can then express several 
cosmologically relevant parameters of the collisional-DM 
model in terms of their collisionless-DM counterparts, such 
as the cosmological redshift, the luminosity distance, and 
the distance modulus of the various light-emitting sources, 
together with the Hubble and the deceleration parameters, 
which characterize the cosmic expansion. 
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3.1. The cosmological redshift 

A supporter of the coUisionless-DM scenario would define 
the corresponding cosmological redshift parameter, z, as 



5 + 1 = 



Rim) 
Riv) ' 



(43) 



which, upon consideration of Eq. (38), is written in the form 

Fivo), 



S + 1 



where 



■(^ + 1) 



.(.o)^fM^l + .... 

Taking into account Ec^. (42), Eq. (44) results in 
1 + wQ M 



S + 1 



l + u>l]M[l + 31n(l + 2)] 



(^ + 1) 



(44) 



(45) 



(46) 



Confining ourselves to relatively low values of the cosmo- 
logical redshift parameter (e.g., z < 5), to ensure that the 
combination Sw^Im hi(l-t-z) remains sufficiently-lower than 
unity even for relatively large values of the combination 
wil.M (e.g., w^Im ^ 0.1, i.e., w ^ -i), we can apply the 
technique used in Eqs. (25) - (27), to obtain 



1 + 5- (1 + z) 



I — SwQm 



(47) 



In this case, we note that, for every (fixed) value of the 
cosmological redshift z, i.e., as defined in the collisional- 
DM model, the corresponding collisionless-DM quantity z 
is always a little bit smaller (z < z). 

In other words, on observing a light-emitting source 
of the coUisional-DM model, an observer who adopts the 
collisionless-DM scenario (realizing redshifts as z instead of 
z), necessarily admits that this source lies a little bit farther 
(z) than expected (z). 

3.2. The luminosity distance and the distance modulus 

Nowadays, the most direct and reliable method for deter- 
mining, observationally, the (relatively) recent history of 
the Universe expansion, is to measure the redshift and the 
apparent luminosity (equivalently, the apparent magnitude, 
to) of cosmologically-distant indicators (standard candles), 
whose absolute luminosity (equivalently, the absolute mag- 
nitude, M) is assumed to be known. 

SN la events constitute one of the most suitable cos- 
mological standard candles (Plionis et al. 2009). With 
the aid of these events, a number of scientific groups 
have attempted to find evidence in support of a recently- 
accelerating stage of the Universe (Garnavich ct al. 1998; 
Riess et al. 1998, 2001, 2004, 2007; Perlmutter et al. 1999; 
Knop et al. 2003; Tonry et al. 2003; Astier et al. 2006; 
Wood-Vasey et al. 2007; Kowalski et al. 2008; Hicken et al. 
2009). In each and every one of these surveys, the SN la 
events, at peak luminosity, appear to be dimmer (i.e., they 
seem to lie farther away) than expected. This result was, 
eventually, accommodated within the context of the con- 
cordance model, by a DE fiuid of negative pressure, with 
nx ^ 0.7 (see, e.g., Spergel et al. 2003; Tonry et al. 2003). 
However, in view of Eq. (47), there may be also another, 
more conventional interpretation. 



Photons travel along null geodesies, ds^ = = ds^, 
which remain unaffected by conformal transformations. 
Accordingly, in both the collisional-DM and collisionless- 
DM approaches, the radial distance of a light-emitting 
source (in comoving coordinates) is the same, i.e.. 



f = c{rir~ rie) = r . 



(48) 



where Tjr and rje are the conformal times of reception and 
emission of light, respectively (usually, r]r = rjo). 

In this case, with the aid of Eq. (47), the formula deter- 
mining the luminosity distance in a spatially-flat collisional- 
DM model 



dL{z)^rSM{l + z) 



(49) 



(see, e.g.. Peacock 1999, p. 92) can be expressed in terms 
of the corresponding collisionless-DM quantity 



Jl(z) =ri?(7?o)(l + 5) 



as 
dL 
dL 



1 



1 + wQ 



M 



-(1 + ^) 



(50) 



(51) 



This relation is very interesting. It suggests that, in a 
Universe containing collisional DM (i.e., as long as w 7^ 0), 
there exists a characteristic (transition) value of the cosmo- 
logical redshift. 



Zc = (1 + w17m)^^^ - 1 



(52) 



such that, the luminosity distance of the various light- 
emitting sources located at z > Zc, is always larger than 
what is inferred by a supporter of the collisionless-DM sce- 
nario. Therefore, an observer who treats the DM as col- 
lisionless (measuring distances in terms of d^) necessarily 
admits that any standard candle located at z > Zc lies far- 
ther than expected, i.e., < dL- 

The same also happens in the case of the distance mod- 
uli corresponding to and d^. The K-corrected distance 
modulus, /i(z) = TO — Af , of a light-emitting source in the 
collisional-DM model is given by 



^l{z) = 5 log 



Mpc 



25 



(53) 



(see, e.g., Narlikar 1983, Eqs. (13.10) and (13.12), p. 359), 
where di is measured in megaparsecs (Mpc). In a similar 
fashion, 



/i(z) = 5 log 



Mpc 



25 



(54) 



is the distance modulus of the same source, as defined by 
someone who, although living in the collisional-DM model, 
insists on adopting the (traditional) coUisionlcss-DM ap- 
proach. Subtracting Eqs. (53) and (54) by parts, and using 
Eq. (51), we obtain 



= fl+ IbwQm log(l + z) — 5 log (1 + wflAi) ■ 



(55) 



According to Eq. (55), any light-emitting source of the 
collisional-DM Universe located at z > Zc, from the point of 
view of an observer who insists on adopting the collisionless- 
DM approach (treating the various distance moduli in 
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terms of /i), appears to be dniimer than expected, i.e., 

We cannot help but notice the prominent similarity be- 
tween the characteristic value Zc and the transition redshift, 
zt , which, according to a supporter of the collisionless-DM 
scenario, signals the onset of the dimming of the SNe la 
standard candles, something that is interpreted (by such an 
observer) as an entry into a phase of accelerated expansion. 
As we see, what actually happens in a coUisional-DM model 
is that, from the point of view of someone who incorrectly 
assumes the DM as collisionless (i.e., measuring cosmolog- 
ical distances in terms of instead of the truly-measured 
quantity di), an inflection point in the versus z diagram 
(namely, Zc) will arise anyway. For w^m =0.1, i.e., w = ^ 
(the DM consists of relativistic particles), the characteristic 
transition value in Eq. (52) is set at Zc = 0.37, while, for 
lower values of w, z^ reaches up to 0.39. These results lie 
within the observationally traced range of values concern- 
ing Zt, namely zt = 0.46 ± 0.13 (Riess et al. 2004, 2007). In 
other words, the spatially-flat collisional-DM model given 
by Eq. (28) does not suffer from the coincidence problem. 

Therefore, if the DM possesses some sort of thermo- 
dynamical content, then, it is possible that (i) the "in- 
fernous" discrepancy between the expected value of the dis- 
tance modulus [Jl) of a SN la standard candle and the cor- 
responding observed one (/j,), and (ii) the accompanying 
inflection point, zt^ that signals the transition from decel- 
eration to acceleration, may both arise only because many 
cosmologists (although living in a collisional-DM model) in- 
sist instead on adopting the (traditional) collisionless-DM 
approach. In the next Section, we demonstrate that this is 
exactly what happens in a collisional-DM Universe. 

3.3. The Hubble and the deceleration parameters 

By virtue of Eq. (38), the Hubble parameter inferred by a 
supporter of the collisionless-DM scenario, H, is written in 
terms of H as 



F S dT]\F 



(56) 



from which, a much more interesting relation can be ob- 
tained in terms of the cosmological redshift. By taking into 
account that 



1 d / 1 \ , , d f 1 , 
Sd-AF)--^'^'^''d-z[F^^ 

Eq. (56) is written in the form 



dz \ F 

which, in view of Eq. (42), results in 
1 — 2wn]\i + Swfli^j ln(l + z) 



H = H 



or else 



(l + «;f7M[l + 31n(l + 2)])2 



(57) 



(58) 



(59) 



H = Fo(l + z)^(i+"'"*') 

1 - 2wQ.M + iwVlM ln(l -I- z) 
^ (l + i«r!M[l + 31n(l + z)])2 ' 



(60) 



where we have also used Eq. (31). We note that, to terms 
linear in w^Mi 



Ho = Hail-4.wnM) , 



(61) 



i.e., within the context of the collisionless-DM approach, 
at the present epoch (when z = 0), the Universe expands 
only as long as wUm < \, and, in any case, at a lower rate 
than the collisional-DM treatment (Hq) implies. In view of 
Eq. (47), i.e., at relatively low values of z, Eq. (60) can be 
written in terms of z, as 



H = iJo (1 + z)2(i--'"'"m) (l-SwriM) 

1 - bw^M + Sw^M ln(l + z) + 0(wf7A/)2 



[1 - 2wnM + SwQm ln(l + z) + OiwilAi)^] 



(62) 



We verify that, to terms linear in wflM, Eq- (61) is (also) 
valid at z = 0. 

By analogy with Eq. (32), a supporter of the 
collisionless-DM scenario would define the corresponding 
deceleration parameter, q, as 



€z^(i+f) 

Hi~z) 



which, by virtue of Eq. (62), yields 



1 



1 - AwQm + Gw^M ln(l + z) + OiwQAif 



1 - lOwf^A/ + 6wnM ln(l + S) + 0(w17a/)2 



(63) 



.(64) 



Now, the condition for accelerated expansion (q < 0) is 
translated to 



[1 - AwflM + GwilM ln(l + z)] X 

[1 - lOwflM + GwriM ln(l + z)]<0, 

from which, to terms linear in w^Im, we obtain 

q{S) < ^ 1 - UiuflM + UwUm ln(l + z) < . 



(65) 



(66) 



From Eq. (66) it becomes evident that, from the point of 
view of someone who insists on adopting the collisionless- 
DM approach, q{z) < at cosmological redshifts 



z < Zt 



1. 



(67) 



This relation is very interesting: It suggests that, if the 
Universe matter-content is treated as a collisional-DM fluid 
with w being larger than a critical value, Wc, such that 



0.0714 



(68) 



(i.e., w > Wc ~ 0.238), then, from the point of view of 
someone who incorrectly assumes the DM as collisionless, 
there exists a transition value, zt, of z, below which, the 
post-recombination Universe (as being realized by such an 
observer) is accelerating, independently of any notion of DE 
or the cosmological constant. 

In other words, if the Universe evolution is driven by a 
collisional-DM fluid with w > Wc, then, the apparent accel- 
eration of the cosmic expansion could (very well) be due to a 
misinterpretation of several cosmologically-relevant param- 
eters, by an observer who (although living in a cosmological 
model filled with collisional DM) insists on adopting the 
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collisionless-DM approach. At the same time, for this ob- 
server, the cosmologicaUy distant indicators would appear 
to be dimmer than expected (cf. Eq. (55)). 

We recall here that the recent observational data con- 
cerning the SNe la standard candles set the transition 
redshift between accelerated and decelerated expansion at 
zt = 0.46 ± 0.13 (Riess et al. 2004). In this case, the com- 
bination of Eqs. (47) and (67) results in the non-linear al- 
gebraic equation involving the transition value, zt, of the 
truly-measured quantity z 



(1 + Zt) eO-25/3«'nA/ ^ 3 2114 (1 + z^f^^-^M 



(69) 



Equation (69) can be solved numerically with respect to the 
combination w^Im- Accordingly, we verify that the value 



{wnM)t = 0.0932 ±0.0060 



(70) 

reproduces (exactly) the above observational result for zt- 
By virtue of Eq. (70), we note that w ~ i, i.e., compati- 
bility of the coUisional-DM approach with the observational 
data, currently available, suggests that, in Eq. (3), the pres- 
sure of the cosmic fluid under consideration is due to radia- 
tion. Nevertheless, in view of Eq. (17), the evolution of the 
rest-mass density indicates that, for every value of w, the 
spatially-flat model given by Eq. (28) is matter-dominated. 
Taken together, these results imply that the DM itself, be- 
ing responsible for the non-vanishing pressure, consists of 
relativistic particles ("hot" DM). We verify this result by 
overplotting Eq. (55) in the Hubble {ji versus z) diagram 
of a SN la dataset. 



3.4. Application to a sample of SN data 

An extended sample of 192 SN la events has been used by 
Davis et al. (2007) to scrutinize the viability of various DE 
scenarios. This sample^ consists of 45 SNe from a nearby 
SN la dataset (Hamuy et al. 1996; Riess et al. 1999; Jha 
et al. 2006), 57 events from SNLS - the Super Nova Legacy 
Survey (Asticr et al. 2006), 60 intermediate- redshift events 
from ESSENCE - the Equation of State: SupErNovae trace 
Cosmic Expansion program (Wood-Vasey et al. 2007), and 
30 high-z SNc from the Gold-07 sample (Riess et al. 2007). 

To ovcrplot Eq. (55) on the ^ versus z diagram of this 
dataset, first of all, we need to determine the functional 
form of the luminosity distance ^^(z) (and, through it, the 
corresponding form of the distance modulus, ii{z)), used 
by someone who, although living in a collisional-DM model, 
insists on the collisionless-DM approach. This observer per- 
forms calculations in the (traditional) framework of a pres- 
sureless Universe, adopting the corresponding formula of 
the luminosity distance. In a spatially-flat model, this for- 
mula is given by (see, e.g., Carroll et al. 1992) 



Ji(z) = ^ (1 + 5)1/^ 



il + ~zf' 



1 



(71) 



representing the luminosity distance in the EdS Universe, 
the (conformally) pressureless counterpart of the 
collisional-DM model given by Eq. (28). 

However, we need to stress that, in a collisional-DM 
Universe, the measured quantity (corresponding to the cos- 
mological redshift) is z and not z, as it is (falsely) admit- 
ted by someone who (mis)trcats the DM as collisionlcss. 



^ Available at http://www.ctio.noao.edu/essence or at 
http : / /braeburn . pha . jhu . edu/ ~ariess/R06 



Therefore, in order to include, also, the function /i(z) in 
the Hubble diagram of the SN la dataset used by Davis et 
al. (2007), we have to express dL{z) in terms of the truly- 
measured quantity, z. It can be done (appropriately) by 
inserting Eqs. (47) and (61) into Eq. (71), to obtain 



2c 



(1 + 2)i(^^^"'^^") - 1 



i(l-3-!uf2Aj) 



(72) 



However, as we have already mentioned, this is not the 
case for a supporter of the collisionless-DM scenario. In de- 
picting Eq. (54) - with ^^(z) being given by Eq. (71) - on 
the fj. versus z diagram of a sample of SN events, this ob- 
server (admitting that w = 0), unavoidably, misinterprets 
the measured quantity z as z (and the quantity Hq as -ffo)- 
In other words, the theoretical formula of the luminosity 
distance that is used by someone who, although living in 
a (spatially-flat) collisional-DM model, insists on adopting 
the collisionless-DM approach is (incorrectly) written in the 
form of Eq. (71) with z simply replaced by z 



(1 



^l/2 1 
z] ' — 1 



(73) 



instead of that given by Eq. (72). In what follows, we admit 
that Hq = 70.5 Km/ sec/ Mpc (Komatsu et al. 2009), hence 
2c/ Ha = 8509.8 Mpc. 

We then ovcrplotted in the Hubble diagram of the SN 
la dataset used by Davis et al. (2007) the theoretical curves 
corresponding to the distance moduli: /i(z) for wflM = 0.16 
(green sohd line); jj.{z) (also for wflM = 0.16) with ^^(z) 
being given by Eq. (72) (orange solid line); and /i(z) with 
dL{z) being given by Eq. (73) (dashed line). The outcome 
is presented in Fig. 1. We observe that, the (appropriately 
translated in terms of z) collisionless-DM quantity jj,{z) - 
with dL{z) being given by Eq. (72) (orange solid line) is 
quite far from being able to reproduce these data, although, 
for z < 1.75, it is much closer to the "real world" (the /i 
versus z distribution of the SN la data available) than the 
incorrectly used quantity jl{z) - with d^i^z) being given by 
Eq. (73) (dashed line). 

The situation changes, completely, when someone takes 
into account the thermodynamical content of a collisional- 
DM fluid with w^M = 0.16, thus using Eq. (55) instead of 
Eq. (54) alone. In this case, the function //(z) (green solid 
line) seems to fit the entire dataset under consideration 
quite accurately. 

As we observe in Fig. 2, apart from a small number of 
SN events (4 of 192), the entire dataset used by Davis et 
al. (2007) lies within the stripe formed by the theoretical 
curves /i(z), corresponding to w^Im = 0.10 (red solid line) 
and wVIm = 0.19 (blue solid line), while the best fit to 
this sample appears to be achieved for w^m = 0.16 (green 
line). In a spatially- fiat cosmological model, it is clear that, 
the collisional-DM treatment is not only much closer to, 
but, actually, can reproduce quite accurately the /i versus 
z distribution of the SN la sample used by Davis et al. 
(2007). 

Nevertheless, there is a "delicate point" in this treat- 
ment. For w^M > 0.10, we have u; > ^. In other words, 
compatibility of the collisional-DM approach with the cur- 
rently available observational data, suggests that the cosmic 
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SN' dala as used bvDaviiel al. (2007) 
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0,0 OJ 10 1.5 2.0 



Fig. 1. Hubble diagram of the SN la sample used by Davis 
et al. (2007). Overplottcd arc the theoretical curves, cor- 
responding to the distance moduli: ii{z) for wD,m — 0.16 
(green solid line); fl{z) (also for wflM = 0.16) with dL{z) 
being given by Eq. (72) (orange solid line); and fi{z) with 
dL{z) being given by Eq. (73) (dashed line). We observe 
that, after the thcrmodynamical content of a coUisional-DM 
fluid is taken into account, the theoretical curve represent- 
ing the distance modulus, fi{z) (Eq. (55)), fits the entire 
dataset quite accurately (green line). 



fluid (i.e., above all, the DM itself) consists of relativistic 
particles. 



SNdataasusedbyDa¥isaal.[200T) 
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Fig. 2. Overplottcd in the Hubble diagram of the SN la 
sample used by Davis et al. (2007), are the theoretical 
curves of the distance modulus in the collisional-DM model, 
i.e., Eq. (55), for several values of the combination wrtM- 
Almost the entire dataset lies within the stripe formed by 
the curves fJ.{z), corresponding to wflM = 0.10 (red solid 
line) and w^Im = 0.19 (blue solid line), while, the best fit to 
this sample is achieved for wi^M = 0.16 (green line). Once 
again, the dashed line represents the (incorrectly used) the- 
oretical curve fl{z) with dL{z) being given by Eq. (73). 



Taken together, these results strongly suggest that, if 
the DM constitutes a perfect fluid (of relativistic particles) 
with thcrmodynamical content, then, what is interpreted 
as the "dimming" of the SNc la standard candles might 
only be apparent, provided that the cosmologists no longer 
insist on adopting the collisionless-DM approach. 

In other words, in a Universe filled with collisional (and 
relativistic) DM, the unexpected dimming of the distant 
light-emitting sources can be explained in a more conven- 
tional way, than that implemented within the context of 
the accelerated expansion. Hence, before inventing any new 
theory, it is useful to allow for a suitable use of the (so 
far) neglected degrees of freedom (energy of the internal 
motions, pressure, etc.). As we have shown, these internal 
physical characteristics can reveal their influence on sev- 
eral parameters of cosmological significance (scale factor, 
cosmological redshift, luminosity distance, Hubble and de- 
celeration parameters) and yield a consistent alternative to 
the currently-accepted DE concept. 

3.5. An independent confirmation from BAO's data analysis 

Among the various (currently available) techniques for trac- 
ing the expansion history of the Universe, baryon acoustic 
oscillations (BAO) appear to have the lowest level of sys- 
tematic uncertainty (Albrecht et al. 2006). 

BAO is a series of peaks and troughs, with wavenumber 
(approximately) O.O6/1 {Mpc)~^ (Eisenstein et al. 2005), 
arising, on large scales, in the power spectrum of matter 
fluctuations after recombination. They occur because the 
primordial cosmological perturbations excite sound waves 
in the relativistic plasma of the early Universe (Silk 1968; 
Peebles & Yu 1970; Sunyaev & Zel'dovich 1970; Bond & 
Efstathiou 1984, 1987; Holtzman 1989; Hu & Sugiyama 
1996). This process continues for a short time-interval af- 
ter recombination, until the epoch known as baryon-drag 
epoch {zd ~ 1089), when, eventually, the baryons are re- 
leased from the Compton "drag" of photons (Eisenstein & 
Hu 1998). 

The BAO's wavenumber is related to the comoving 
distance of the sound horizon at the baryon-drag epoch, 
Tg, which depends on both the total- mass density of the 
Universe and the corresponding baryonic quantity (see, 
e.g., Wang 2006). WMAP constraints on n^jh'^ and fle/i^ 
(Komatsu et al. 2009) suggest that rs{zd) — 153.5 Mpc. In 
other words, BAO are met on relatively-large scales, which, 
at the present epoch, are (still) in the linear regime. It is, 
therefore, expected that, these acoustic signatures should 
be present (also) in the distribution of galaxies (Goldberg 
& Strauss 1998; Meiksin, White & Peacock 1999; Seo & 
Eisenstein 2005; Eisenstein, Seo & White 2007). 

Since the comoving distance of the sound horizon at the 
baryon-drag epoch is known, we can use BAO as standard 
rulers, to determine the functional form of the distance - 
redshift relation in the Universe. In particular, by determin- 
ing observationally the apparent size of a BAO's peak (it is 
identified as a clustering within a galaxy distribution), we 
can extract reliable estimates of the Hubble parameter, as 
well as the angular diameter distance, at the cosmological 
redshift, in which, this acoustic signature is observed (Blake 
& Glazebrook 2003; Hu & Haiman 2003; Sco & Eisenstein 
2003). 
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The spectroscopic Sloan Digital Sky Survey (SDSS) 
Data Release 7 (DR7) represents the final set of galax- 
ies, observed using the original SDSS-target selection cri- 
teria (York et al. 2000). When combined with the data re- 
leased from the 2-degree Field (2dF) galaxy-redshift survey 
(Colless et al. 2003), the resulting sample comprises 893319 
galaxies and covers a solid angle of 9100 deg^. An analy- 
sis of the clustering of galaxies within this sample, inferred 
BAO signals in the power spectrum measured in several 
slices of the cosmological redshift (Percival et al. 2010). 

To minimize the systematic errors arising from treat- 
ing the linc-of-sight dilation in an equivalent manner to 
the transverse one (a common bias, involving the so-called 
Alcock-Paczyhski (1979) effect), cosmologists defined the 
distance measure 



Dv{Z) 



II' 1 



{l + ZfD\{Z) 



cZ 

Wz) 



1/3 



(74) 



(Eisenstein et al. 2005; Percival et al. 2007), where Z 
is the observationally-determined (i.e., spectroscopically- 
measured) value of the cosmological redshift of a large-scale 
structure, and 



Da{Z) = 



Sim) 
i + z 



(75) 



is the corresponding angular diameter distance in a 
spatially-fiat cosmological model (see, e.g.. Peacock 1999, 
p. 22). The combination of Eqs. (74) and (75) results in 



DviZ) 



, cZ 

'h{z) 



1/3 



(76) 



where, once again, r is the radial distance (in comoving co- 
ordinates) of the structure formation under consideration. 

According to Percival ct al. (2010), there exists a robust, 
statistically-independent (distance) constraint, arising from 
BAO's data, which involves the value of the ratio 



/ 



i:'v(0.35) 
Dv{0.2) 



1.736 ±0.065. 



(77) 



Clearly, in a collisional-DM model, Eq. (77) represents the 
truly-measured value of /. In this case, we ask ourselves 
what is observed by someone who, although living in a 
collisional-DM model, insists on adopting the (traditional) 
collisionless-DM approach. 

To answer this question, we note that, when, cither 
the coUisional-DM-oricnted observer or the corresponding 
collisionless-DM one, spectroscopically measures the cos- 
mological redshift of a particular large-scale structure, they 
both refer to the same quantity, Z. Accordingly, a supporter 
of the collisionlcss-DM scenario would express Eq. (76) in 
the form 



DviZ) 



, cZ 



1/3 



(78) 



where we have also used Eq. (48). Equation (78) differs from 
Eq. (76) only in the definition of the present-time value of 
the scale factor and in the functional dependence of the 
Hubble parameter on Z. In this case, upon consideration of 



Eqs. (45) and (60), we obtain 



DvjZ) 
DviZ) 

DvjZ) 
DviZ) 



W^m)'^' ( ^ 



1/3 



(l + wf^A/) 



IT 

2/3 [1 -f- wVLm 



(79) 

■iw^M ln(l + Z)f'^ 



[1 - 2wnM + '^wUm ln(l + Z)] 



1/3- 



Now. as far as a supporter of the collisionlcss-DM scenario 
is concerned, the BAO's constraint given by Eq. (77) is 
translated as 



£'y(0.35) 
bviQ.2) 



1 - 2wnM + iwVlM ln(1.2) 1 ^'^ 
1 - 2wnM + iw^M ln(1.35) 

l + wnM + SwQm ln(1.35)l 



which for wfl 



M 



1 -I- wQ M 
DviQ.35) 
DviO.2) ' 

0.10 results in 



iwilM ln(1.2) 



i)y(0.35) 



= 1.766. 



^~^'"''"="-' DviO.2) 
and for wCIm = 0.16, we have that 



/I 



w^1m—0. 16 



Dy iO.35) 
DviO.2) 



1.788. 



(80) 



(81) 



(82) 



Both values of / lie within the range of values [1.671 , 1.801] 
of the measured quantity, /, given by Eq. (77). Therefore, 
as far as a supporter of the collisionless-DM scenario is 
concerned, the observational constraint, /, inferred from 
BAO's data, is (at least) compatible with (if not almost 
identical to) the corresponding constraint, /, obtained 
within the context of the collisional-DM model. 

However, when the coUisionless-DM-oricnted observer 
attempts to verify Eq. (81) (or Eq. (82)) theoretically, a 
controversy arises. This observer applies Eq. (78) to the EdS 
Universe (the pressureless counterpart of the collisional- 
DM model under study), which accommodates the stan- 
dard iw = 0) cold dark matter iSCDM) cosmology. In 
this case. 



DviZ) = 2^/^^Z^^^^ — — , 

^ ' Hq 1 + Z5/6 



(83) 



hence, a supporter of the collisionless-DM scenario obtains 



fscDM 



£'v(0.35) 
bviQ.2) 



1.553 



(84) 



(see, also, Percival et al. 2010). Clearly, there is a difference 
between the theoretical result of Eq. (84) and the obser- 
vational one given by Eq. (81) (or Eq. (82)). As far as the 
supporter of the collisionless-DM scenario is concerned, one 
way to compensate this difference, is to impose the exis- 
tence of an extra (dark) energy component or the cosmo- 
logical constant. However, even when he/she does so, i.e., 
within the context of the KCDM model, the theoretical 
value of the distance constraint, /, induced by the signa- 
ture of BAO on cosmic structure, reaches up to 



fhCDM = 1-670 



(85) 
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(Eisenstein et al 2005; Percival et al. 2010). Although 
marginally, /acdm lies outside the range of values of the 
truly-measured quantity, /, given by Eq. (77). In view of 
the above-mentioned results, we may conclude that, within 
the context of the coUisional-DM approach, the BAO's- 
oriented observational data require (and acquire) a more 
sophisticated interpretation than the one provided by the 
(collisionless) KCDAI model. 

4. Discussion 

We have examined the possibility that the extra (dark) en- 
ergy needed to flatten the Universe is represented by the 
energy of the internal motions of a collisional-DM fluid. 
Accordingly, we have considered the evolution of a cosmo- 
logical (toy-) model driven by a gravitating fluid (consisting 
of DM - dominant - and baryonic matter - subdominant) 
with thermodynamical content. As a consequence, the en- 
ergy of this fluid's internal motions has also been taken 
into account as a source of the universal gravitational field. 
Accordingly, we have asked ourselves, whether this model 
can also accommodate the apparent dimming of the cos- 
mologically distant indicators and the associated phase of 
accelerated expansion. 

In particular, since observational data indicates that 
WIMPs (of which the DM is believed to consist) can 
be collisional, we have assumed that the matter of the 
Universe (although resembling test particles receding from 
each other) can be represented by the volume elements of 
a (classical) collisional-DM fluid with some sort of inter- 
nal structure, hence thermodynamical content. In this way, 
we have been able to determine the "correct" form of the 
scale factor, which (under the assumption that the DM is 
thermodynamically involved) governs the evolution of the 
Universe (modeled as a spatially-flat RW space-time), in 
addition to a series of parameters of cosmological signifi- 
cance. 

Our findings are quite promising. In principle, the en- 
ergy of the internal motions of the collisional-DM fluid can 
account for the (extra) DE, so that, at the present epoch, 
J7 = 1 (cf. Eq. (23)), while the post- recombination Universe 
remains ever-decelerated (cf. Eq. (33)). 

We next attempted to determine what is inferred by 
someone who, although living in a collisional-DM model, 
insists on adopting the (traditional) collisionlcss-DM ap- 
proach. 

To do so, we have applied the technique developed by 
Kleidis & Spyrou (2000). With the aid of this technique, 
we have derived the (conformal) transformation (cf. Eqs. 
(38) and (42)), which relates the collisional-DM descrip- 
tion of a cosmological model (in terms of which p, 11 > 
and ^ 7^ 0) to the corresponding coUisionless-DM (pres- 
sureless) approach. With such a "tool" at hand, we have 
explored the way that a supporter of the collisionless-DM 
scenario interprets observations carried out in a collisional- 
DM Universe. In passing, we note that, within the context 
of general relativity, every effort to treat a collisional-DM 
model as pressureless is questioned (Appendix B). 

The debate between collisional- and collisionless-DM 
approach is, definitely, in favor of the former. In partic- 
ular, for every value of the cosmological redshift (z), as it 
is defined in the collisional-DM model, the corresponding 
collisionless-DM quantity, z, is always a little bit smaller (cf. 



Eq. (47)). As a consequence, in the collisional-DM model 
there is a characteristic value of the cosmological redshift, 
Zc (cf. Eq. (52)), above which, the luminosity distance of the 
various light-emitting sources is always higher than what is 
inferred by an observer who treats the DM as pressureless 
(cf. Eq. (51)). In other words, from the point of view of 
someone who (although living in a collisional-DM model) 
insists on adopting the (traditional) collisionless-DM ap- 
proach, the cosmologically-distant indicators, located at 
z > Zci seem to lie farther away (i.e., appear dimmer) than 
expected (cf., also, Eq. (55)). 

The similarity between the characteristic value z^ and 
the (observationally-traccd) transition redshift, Zt, which, 
according to a supporter of the collisionless-DM scenario, 
signals the onset of the dimming of the SNe la standard 
candles, is obvious. 

On the other hand, after the thermodynamical con- 
tent of a collisional-DM fluid is taken into account, the 
theoretical curve representing the distance modulus, [i(z) 
(now given by Eq. (55)), fits the Hubble diagram of an ex- 
tended sample of SN la standard candles quite accurately 
(green solid line in Fig. 1), in contrast to the correspond- 
ing collisionless-DM quantity, \^.(z)^ given either (appropri- 
ately) by the combination of Eqs. (54) and (72) (orange 
solid line in Fig. 1) or (incorrectly) by the combination of 
Eqs. (54) and (73) (dashed line in Fig. 1). In fact, a cosmo- 
logical model filled with collisional DM could accommodate 
the majority of the currently-available observational data, 
including, also, those related to BAO, without the need for 
any dark energy or the cosmological constant (cf. Eqs. (81) 
and/or (82)). 

At the same time, from the point of view of an observer 
who treats the DM as collisionless, the Universe appears 
to be either accelerating or decelerating, depending on the 
value of the cosmological redshift (cf. Eq. (67)). 

In this case, the quantity w, which, in the collisional- 
DM approach, parameterizes the various flows, also plays 
another (more interesting) role. As we have found, for 
w^M > 0.0714, there exists a (theoretically-determined) 
transition value, zt, of the (collisionless-DM-oriented) cos- 
mological redshift, z, such that, for z < Zj, we have g < 0, 
i.e., from the point of view of someone who adopts the (tra- 
ditional) collisionless-DM approach, the Universe is acceler- 
ating, without the need for any DE or the cosmological con- 
stant. Accordingly, taking into account the observational re- 
sult that the transition redshift between accelerated and de- 
celerated expansion is set at the value Zt = 0.46±0.13 of the 
truly measured quantity z, we have determined the precise 
value of the combination wJIa/, for which the collisional- 
DM approach to the post-recombination Universe is com- 
patible with observations, i.e., (w^M^t ~ 0-0932 ± 0.0060. 
This result is in complete agreement with what is implied 
by Fig. 2, i.e., it; ~ i. In other words, compatibility of 
the collisional-DM treatment under study with the observa- 
tional data currently available, suggests that the DM itself 
consists of relativistic particles (hot DM). 

In conclusion, the assumption that the DM constituents 
can be both collisional and relativistic, could provide a 
reasonable and conventional explanation for several open 
aspects of modern cosmology, including: (i) The extra 
(dark) energy needed to flatten the Universe, which can 
be compensated by the energy of the internal motions of 
the collisional-DM fluid, (ii) The observed dimming of 
the SNe la standard candles and the apparent acceler- 
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ated expansion of the Universe, both of which might be 
due to the misinterpretation of several cosmologically rel- 
evant parameters by those observers who, although liv- 
ing in a collisional-DM Universe, insist on adopting the 
collisionless-DM approach. The absence of any ad hoc fine- 
tuning in our study, makes the above (purely theoretical) 
results very promising. 

In spite of all the above advantages of our model, we 
have to point out that this is, definitely, a toy model. The 
collisional-DM approach, developed in this article, can (and 
should) be debated, at least, along the following lines: 

(i) The collisional-DM treatment of the Universe does 
not alleviate the age problem^ but, rather, makes it harder. 
In the model given by Eq. (28), the coordinate time, t, is 
related to the corresponding conformal quantity (77) by 



t 



S{r])d7] 



1 



/ 2 



1 + w^Im \iHo 



3(i+mnjtj) 



(86) 



Admitting that, at the present epoch (when t = to and 
rj = rio) Ho w 70.5 Km/sec/Mpc (Komatsu ct al. 2009), 
we find that, within the context of the model given by Eq. 
(28), the age of the Universe is 



in = 



1 



1 + wQ 



M 



X 9.3 Gys , 



(87) 



i.e., less than 9.3 billion years. Clearly, Eq. (87) could be a 
serious drawback of the spatially-flat collisional-DM model. 

(ii) The compatibility of the collisional-DM approach 
with the observational data currently available (cf. Eq. (70) 
and/or Fig. 2), suggests that the matter content of the dark 
sector consists of hot DM. For the time being, the conven- 
tional theory of hot DM does not appear to conform with 
the large-scale structure of the Universe (see, e.g.. Hooper 
2009), although there are recent studies that appear to chal- 
lenge this result (see, e.g., Farrar & Peebles 2004; Gubser & 
Peebles 2004). Nevertheless, a hot-DM model looks much 
less exotic than most of the (currently-investigated) DE sce- 
narios. 

In any case, the assumption that the Universe matter 
content (basically its DM component) can be collisional 
(in the sense that it also possesses some sort of thermody- 
namical content), is to be seen as a natural effort to take 
into account all the (so far, practically, neglected) internal 
physical characteristics of a classical cosmological fluid as 
sources of the universal gravitational field. 

As we have shown, under this assumption, one can com- 
pensate for the majority of the recent observational data, 
inferring that $7 ~ 1, as well as the "unexpected" dimming 
of the SNe la standard candles and the " apparent" acceler- 
ated expansion of the Universe, not to mention the compat- 
ibility with BAO's data analysis, without the need for any 
DE or the cosmological constant, and (certainly) without 
suffering from the coincidence problem. From the point of 
view of someone who (although living in the collisional-DM 
model) insists on adopting the collisionless-DM approach, 
an inflection point (in the versus z diagram) arises, any- 
way, at relatively low values of the cosmological redshift. 

Although speculative, the idea that the DE (needed to 
flatten the Universe) could be attributed to the internal 
motions of a collisional-DM fluid, is (at least) intriguing and 
should be explored further and scrutinized in the search for 
conventional alternatives to the DE concept. 
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Appendix A 

In this Appendix, we present the exact solution to Eq. (24) 
in the text and its subsequent reduction to Eq. (28). In 
paricular, upon consideration of Eq. (12), Eq. (24) reads 



1 + 3wQ. M In 



and can be cast in the (more convenient) form 
2 




HqSq 



2 r 



1 + SwflM In 



(A2) 



In an expanding Universe, Eq. (A2) is valid for every 77 
within the past light-cone (where S < Sq) and it can be 
solved in terms of the error function (see, e.g., Abramowitz 
& Stegun 1970, p. 297) as 



\/6wQ,M 



^l/6'ujf2j\/ 



6wn 



M 



Equation (A3) becomes particularly transparent (and use- 
ful) if we restrict ourselves to the limiting case, where 

For wQm <SC 1, the argument of the error function in Eq. 
(A3) becomes very large. In this case, the special function 
under consideration behaves as 



erf{ip) 



1 



1 1 



(^4) 



(see, e.g., Lebedev, Eq. (2.2.4), p. 19) and, therefore, Eq. 
(A3) results in 



(4) 



1 + 3i«17A/ln(|i) 



(1 + 3wO 
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We now take the natural logarithm on both sides of Eq. 
(A5), to obtain 
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Within the post-recombination era, ^ < 1090, hence 

In(f^)^ < 21. Therefore, as long as wf^M < 1, Eq. (26) 
holds, hence Eq. (A6) is written in the form 



which, in view of the definition given by Eq. (29), results 
in 

2 

l + 3mS2,, 



that is, Eq. (28). 



Appendix B 



In curved space-time, the Einstein-Hilbert (EH) action, 
which governs the dynamical evolution of the gravitational 
field, 5^1/, within the context of GR, is given by (see, e.g., 
Papapetrou 1973, Eq. (33.29), p. 122) 



(Bl) 



where g is the determinant of the metric tensor g^^u, 
is the scalar curvature, k ~ '^■nGjc^^ and T is the trace 
of the energy-momentum tensor of the Universe matter- 
energy content. In a similar fashion, the action of the grav- 
itational field in terms of 5^^, i.e., after the transformation 
indicated in Eq. (36) is performed, is written in the form 



g 7^ + kT d a; 



which, upon consideration of Eqs. (36) and (39) of Kleidis 
& Spyrou (2000), results in 



AF 



'-EH 



(S3) 



where AF = g^'^F.^ is the d' Alembert operator with re- 
spect to the original metric, g^v It is evident that, every 
effort to express X in terms of the original gravitational 
field, yields the appearence of extra terms in the grav- 
itational action, in addition to the EH one. In this way, 
variation of Eq. (B3) with respect to g^'^ does not lead to 
the Einstein field equations of GR (i.e., those derived from 
Eq. (Bl)), thus (necessarily) resulting in a modified theory 
of gravity (in connection, see, e.g., Bruneton & Esposito- 
Farese 2007). 
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